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I . Introduction 

This report contains a target position estimation program for 
each of the following calculators: The Hewlett-Packard HP-41CV, 

the Sharp PC-1500 (or Radio Shack TRS-80 PC-2) , the Sharp PC-1211 
(or Radio Shack TRS-80 PC-1) , the Casio FX-702P and the Texas 
Instruments TI-59. 

The programs provide a means of implementing a position esti- 
mation procedure that is described in Appendix 1. The procedure 
is based on the following assumptions: Bearings taken on or from 

a target are available from two or more stations of known position. 
The positions of the stations and the target are such that they 
can be considered to be on the surface of a plane (a flat earth 
model) . The error in a bearing taken on or from a station is a 
normal random variable. Its standard deviation (bearing error) 
is known and its mean (bias) is zero (if there is bias in a bear- 
ing, it is known and it is removed). The errors in bearing mea- 
surements are independent. 

The position estimation procedure requires an initial estimate 
of the target's position. In the programs that are given here, 
the initial estimate is at the intersection of the bearing lines 
determined by the first two target bearings that are input to 
the program. Because of the use of this method to choose an ini- 
tial estimate, the first two stations, with respect to the order 
of data input, should be chosen so that the intersection of their 
bearing lines is likely to be closer to the target's position 
than the intersection of the bearing lines from any other pair 
of stations. 
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The target ranges and errors of the bearings of the first two 
stations determine the distance between the intersection of their 
bearing lines and the target's position. In particular, if the 
angular separation between two stations as seen from the target 
is small relative to the bearing error of one or both of the sta- 
tions, the bearing lines from the two stations may not intersect. 
If they do not intersect and they are not parallel, their reci- 
procal bearing lines will. In this case, if the observed bearings 
from the two stations were the first two bearing inputs to the 
program, the initial estimate of the target's position would be 
at the intersection of the two reciprocal bearing lines, and a 
gross error in the final position estimate could result. 
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II. Program User Instructions 

Before using a program for the first time, refer to the notes 
that follow the user's instructions. Also, note the comments in 
Section I regarding the relationship between the accuracy of the 
computed estimates and the order of data entry. 

For a PC-1211, use the PC-1500 User Instructions. Use the 
DEF mode and substitute SHFT for DEF wherever DEF appears in the 
instructions . 

Station positions are determined with respect to a known ref- 
erence point. The reference point can be a station position, in 
this case the station bearing and range from the reference point 
are both zero. 

Program listings are given in Section IV and two example 
applications are discussed in Section III. 
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Position Estimation Program 
HP-41CV User Instructions 



Step 


Instruction 


Prompt 


Press 


1 


Select USER mode and enter the program 
(see Note 1). Press £+ to run the 
program. 




E + 


2 


Key in, in decimal degrees, the obser- 
ved bearing of the target from a 
station or the reciprocal of the ob- 
served bearing of the station from 
the target (see Note 2) . 


OBS BRG? 


R/S 


3 


Key in the station bearing in decimal 
degrees from the reference point. Use 
zero if the station is the reference 
point. 


STA BRG? 


R/S 


4 


Key in the station range (in any units) 
from the reference point. Use zero 
if the station is the reference point. 


STA RNG? 


R/S 


5 


Key in the bearing error (standard 
deviation) in decimal degrees. 




R/S 


6 


Repeat Steps 2, 3, 4 & 5 for one or 
more additional stations. 






7 


Compute bearing & range estimates. 




1/x 


8 


A computed target bearing estimate in 
decimal degrees is displayed. 


BRG=e s t • 


R/S 


9 


A computed target range estimate in 
meters. For an elliptical contain- 
ment region, go to Step 10 or Step 17. 
To enter additional data from new or 
old stations, go to Step 24. 


RNG=est . 




10 


For a containment ellipse of a given 
containment probability, press LN. 




LN 


11 


Key in the desired containment proba- 
bility. 


PRB? 


R/S 


12 


Computed value of the ellipse size: 
(See Note 3.) 


SIZE=val . 


R/S 
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Step 


Instruction 


Prompt 


Press 


13 


Computed semi-major axis length: 


SMJ=val . 


R/S 


14 


Computed major axis direction: 


DRC=val . 


R/S 


15 


Computed semi-minor axis length: 


SMI=val . 


R/S 


16 


Computed containment ellipse area: 
(See Note 4 . ) 


A = val. 




17 


For a containment ellipse of a given 
size, press LOG. 




LOG 


18 


Key in the desired containment ellipse 
size (see Note 3) . 


SIZE? 


R/S 


19 


Computed value of the containment 
probability: 


PRB=val . 


R/S 


20 


Computed semi-major axis length: 


SMJ=val . 


R/S 


21 


Computed major axis direction: 


DRC=val . 


R/S 


22 


Computed semi-minor axis length: 


SMI=val. 


R/S 


23 


Computed containment ellipse area: 
(See Note 4.) 


A = val. 




24 


To enter additional data from new or 
old stations, press /x and then repeat 
Steps 2, 3, 4 & 5 
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Notes : 



1. The program size is 35. The key assignments are: TPE -*■ £+, 

CON -*• /x, EST -*■ 1/x, SIZ -»• LOG and PRB LN. If they are not 
present, they must be made in order to follow the user 
instructions. For an alternative, see Appendix 2. 

2. Reciprocal bearings are used when bearings are taken on known 
positions from the unknown position (target) . 

3. In the model that is the basis for the program, for a given 
probability of containment, the minimum area containment region 
is an ellipse centered on the position estimate. The semi- 
major axis = ka„._ and semi-minor axis = ka.._ where k is the 

size of the ellipse and c,, T and a„_ are the standard deviations 

MJ MI 

(uncertainty measure) of the position estimate in the major 
axis and minor axis directions. 

4. The area units are the range units squared. 
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Position Estimation Program 
PC-1500 User Instructions 



; .tep 


Instruction 


Prompt 


Press 


1 


Enter the program. To run the program, 
press DEF, A. For a PC-1211, see Note 5. 




DEF A 


2 


Key in, in decimal degrees, the observed 
bearing of the target from a station or 
the reciprocal of the observed bearing 
of the station from the target (see 
Note 1) . 


OBS BRG? 


ENTER 


3 


Key in the station bearing in decimal 
degrees from the reference point. Use 
zero if the station is the reference 
point. 


STA BRG? 


ENTER 


4 


Key in the station range (in any units) 
from the reference point. Use zero if 
the station is the reference point. 


STA RNG? 


ENTER 


5 


Key in the bearing error (standard 
deviation) in decimal degrees. 


BRG ERR? 


ENTER 


6 


Repeat Steps 2, 3, 4 & 5 for one or 
more additional stations. 






7 


Compute bearing & range estimates. 




DEF Z 


8 


A computed target bearing estimate in 
decimal degrees is displayed. 


BRG=est . 


ENTER 


9 


A computed target range estimate is dis- 
played. For an elliptical contain- 
ment region, go to Step 10 or Step 17. 

To enter additional data from new or 
old stations, go to Step 24. 


RNG=est . 




10 


For a containment ellipse of a given 
containment probability, press DEF, X. 




DEF X 


11 


Key in the desired containment proba- 
bility. 


PRB? 


ENTER 


12 


Computed value of the ellipse size: 
See Note 2. 


SI ZE=val . 


ENTER 
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Step 


Instruction 


Prompt 


Press 


13 


Computed semi-major axis length: 


SMJ=val . 


ENTER 


14 


Computed major axis direction: 


DRC=val . 


ENTER 


15 


Computed semi-minor axis length: 


SMI=val . 


ENTER 


16 


Computed containment ellipse area: 
(See Note 3 . ) 


A = val. 




17 


For a containment ellipse of a given 
size, press DEF, S. 




DEF S 


18 


Key in the desired containment ellipse 
size (see Note 2) . 


SIZE? 


ENTER 


19 


Computed value of the containment 
probability : 


PRB=val . 


ENTER 


20 


Computed semi-major axis length: 


SMJ=val . 


ENTER 


21 


Computed major axis direction: 


DRC=val . 


ENTER 


22 


Computed semi-minor axis length: 


SMI = val . 


ENTER 


23 


Computed containment ellipse area: 
(See Note 3 . ) 


A = val. 




24 


To enter additional data from new or 
old stations, press DEF, C and then 




DEF C 




repeat Steps 2, 3, 4 & 5. 




u 
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Notes : 



1. Reciprocal bearings are used when bearings are taken on known 
positions from the unknown position (target). 

2. In the model that is the basis for the program, for a given 
probability of containment, the minimum area containment 
region is an ellipse centered on the position estimate. The 
semi-major axis = ka_._ and semi-minor axis = kcr _ where k is 
the size of the ellipse and o. An . and a.. T are the standard 
deviations (uncertainty measure) of the position estimate 

in the major axis and minor axis directions. 

3. The area units are the range units squared. 

4. For a definition of the program initiating keys and their 
function, press DEF, H. The display will show: 

TPE = A EST = Z SIZ = S PRB = X. Next press ENTER. The 
display will show: CON = C. To repeat the display, press 
ENTER. 

5. For a PC-1211, use the DEF mode and substitute SHFT for 
DEF wherever DEF appears in the instructions. 
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Position Estimation Program 



FX-702P User Instructions 



S 



Step 


Instruction 


Prompt 


Press 


1 


Enter the program (see Note 1) . To run 
the program, first press FI, 0 . 




FI 0 


2 


Next, key in 1 and press EXE. 


OPTION? 


1 EXE 


3 


Key in, in decimal degress, the obser- 
ved bearing of the target from a 
station or the reciprocal of the ob- 
served bearing of the station from 
the target (see Note 2) . 


OBS BRG? 


EXE 


4 


Key in the station bearing in decimal 
degrees from the reference point. Use 
zero if the station is the reference 
point . 


STA BRG? 


EXE 


5 


Key in the station range (in any units) 
from the reference point. Use zero 
if the station is the reference point. 


STA RNG? 


EXE 


6 


Key in the bearing error (standard 
deviation) in decimal degrees. 


BRG ERR? 


EXE 


7 


Repeat Steps 3 , 4 , 5 & 6 for one or 
more additional stations. 






8 


To compute bearing & range estimates, 
first press FI, 0 . 




FI 0 


9 


Next, key in 2 and press EXE. 


OPTION? 


2 EXE 


10 


A computed target bearing estimate in 
decimal degrees is displayed. 


BRG— est . 


CONT 


11 


A computed target range estimate is dis- 
played. For an elliptical contain- 
ment region, go to Step 12 or Step 20. 

To enter additional data from new or 
old Stations, go to Step 28. 


RNG=est . 




12 


For a containment ellipse of a given 
containment probability, first press 
FI, 0 . 




FI 0 


13 


Next, key in 3 and press EXE. 


OPTION? 


3 EXE 
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Instruction 



Prompt 



Key in the desired probability. 

Computed value of the ellipse size: 

See Note 3. 

Computed semi-major axis length: 

Computed major axis direction: 

Computed semi-minor axis length: 

Computed containment ellipse area: 

See Note 4. 

For a containment ellipse of a given 
size, first press FI, 0 . 

Next, key in 4 and press EXE. 

Key in the desired containment ellipse 
size (see Note 3) . 

Computed value of the containment 
probability: 

Computed semi-major axis length: 
Computed major axis direction: 

Computed semi-minor axis length: 

Computed containment ellipse area: 

See Note 4. 

% 

To enter additional data from new or 
old stations, first press FI, 0 . 

Next, key in 5, press EXE and then 
repeat Steps 3, 4, 5 & 6. 



PRB? 

SIZE=val . 

SMJ=val . 
DRC-val . 
SMI=val . 

A = val. 



OPTION? 

SIZE? 

PRB=val . 

SMJ=val . 
DRC=val . 
SMI=val . 
A - val. 



OPTION? 



Press 

EXE 

CONT 

CONT 

CONT 

CONT 

FI 0 

4 EXE 
EXE 

CONT 

CONT 

CONT 

CONT 

FI 0 

5 EXE 
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Notes : 



1. Enter the program in P0 for FI 0 activation. Before running 
the program, first press F2 then DEFM then 1 and then EXE. 
This is required in order to use the array variables in the 
program. 

2. Reciprocal bearings are used when bearings are taken on known 
positions from the unknown position (target) . 

3. In the model that is the basis for the program, for a given 

probability of containment, the minimum area containment 

region is an ellipse centered on the position estimate. The 

semi-major axis = ka,„_ and semi-minor axis = ka„_ where k 

is the size of the ellipse and and o., T are the standard 

r MJ MI 

deviations (uncertainty measure) of the position estimate 
in the major axis and minor axis directions. 

4. The area units are the range units squared. 

5. For a definition of the program options and their function, 

press FI, 0. Then after OPTION? is displayed, enter 0 and 
press EXE. The display will show: TPE = 1 EST = 2 SIZ = 3. 

Next, press CONT. The display will show: PRB = 4 CON = 5. 

To repeat the displays, press CONT. 
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POSITION ESTIMATION PROGRAM 



TI-59 USER INSTRUCTIONS 



— ■■ 

Step 


Instruction 


Enter 


Press 


Display 


1 


Enter the program (see Note 1) . 








2 


To run the program, press A. 




A 


9 


3 


Enter, in decimal degrees, the observed 
bearing of the target from a station or 
the reciprocal of the observed bearing 
of the station from the target (see 
Note 2) . 


0 


B 


0 


4 


Enter the station bearing in decimal 
degrees from the reference point. Use 
zero if the station is the reference 
point . 


a 


R/S 


a 


5 


Enter the station range from the refer- 
ence point. Use zero if the station is 
the reference point. Use any units. 


P 


R/S 


P 


6 


Enter the bearing error (standard devi- 
ation) in decimal degrees. After press- 
ing R/S , the display indicates the order 
number of the data entry. 


e 


R/S 


n 


7 


Repeat Steps 3, 4, 5 & 6 for at least 
one more station. 
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A 

Display a bearing estimate <p with respect 
to the reference station. 




C 


/\ 


9 


Display a range estimate r with respect 
to the reference station. For an ellip- 
tical containment region, go to Step 10 
or Step 15. To enter additional data 
from new or old stations, repeat Steps 3, 
4, 5 & 6. 




R/S 


/N 

r 


.0 


For a containment ellipse of a given con- 
tainment probability, enter the contain- 
ment probability. Next, press E and 
display the ellipse size (see Note 3) . 


P 


E 


k 


.1 


Display the semi-major axis length. 




R/S 


ko MT 

MJ 


.2 


Display the major axis direction. 




R/S 


Y 


.3 


Display the semi-minor axis length. 




R/S 


ko,„ T 

MI 
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Step 


Instruction 


Enter 


Press 


14 


Display the containment ellipse area. 
(See Note 4.) 




R/S 


15 


For a containment ellipse of a given 
size, enter the containment ellipse 
size (See Note 3) . Next, press D and 
display the containment probability. 


k 


D 


16 


Display the semi-major axis length. 




R/S 


17 


Display the major axis direction. 




R/S 


18 


Display the semi-minor axis length. 




R/S 


19 


Display the containment ellipse area. 
(See Note 4.) 




R/S 


20 


To enter additional data from new or 
old stations, repeat Steps 3, 4, 5 
& 6. 







Displ 



area 



P 



ka 

MJ 

Y 

ka.._ 

MI 

area 
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Notes : 



1. The program requires the normal partition. If the calculator 
has been in use, this can be assured by turning the calculator 
off and then on before loading the program. 

2. Reciprocal bearings are used when bearings are taken on known 
positions from the unknown position (target) . 

3. In the model that is the basis for the program, for a given 
probability of containment, the minimum area containment region 
is an ellipse centered on the position estimate. The semi- 
major axis = ka,, T and semi-minor axis = ka.. T where k is the 

, J MJ MI 

size of the ellipse and a,,, and are the standard deviations 

r MJ MI 

(uncertainty measure) of the position estimate in the major 
axis and minor axis directions. 

4. The area units are the range units squared. 

5. Negative bearing estimates and negative major axis directions 
can result. To convert a negative bearing estimate to the 
value that would be output by the other programs, add 360° to 
the estimate. For example, -5° becomes 355°. To convert a 
negative direction, add 180°. For example, -5° becomes 175°. 
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III. Two Examples 

In Scenario 1, the scenario for the first example, bearings 
are taken on a target from three separate stations. Figure 1 on 
Page 18 shows Scenario 1 and Table 1 below gives the station data. 
The stations are numbered according to the order of station data 
input to the program. 







TABLE 1 








OBS BRG 


STA BRG 


STA RNG BRG 


ERR 


Station 1 


038° 


334° 


13500 


4° 


Station 2 


324° 


050° 


11350 


3° 


Station 3 


003° 


000° 


0 


4° 


Note from Table 


1 that the 


reference 


point is at Station 3 





Program outputs for Scenario 1 are given in List 1 on Page 
19. List 1 and List 2 which gives program outputs for Scenario 
2, are copies of printer tapes that were generated with a Casio 
FP-10 printer and a Casio FX-702P using the Casio FX-702P program. 
With allowance for differences in user instructions, display for- 
mats and round-off errors, the tapes indicate the output that 
should be obtained using any of the other calculator programs ex- 
cept for the TI-59 program where equivalent negative angles occur. 

In Scenario 2, the scenario for the second example, bearings 
are taken from a target on three separate stations. Figure 2 on 
Page 20 shows Scenario 2 and Table 2 below gives the station data. 
As in Scenario 1, the stations are numbered according to the order 
of station data input to the program. In this scenario the reci- 
procal of the observed bearings are used in order to provide an 
equvalent scenario that is appropriate for the program. 
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TABLE 2 



Station 1 
Station 2 
Station 3 



OBS BRG 
211 ° 
172° 
146° 



RCP BRG 
031° 
352° 
326° 



STA BRG 
000 ° 
115° 
082° 



STA RNG 
0 

11100 

13800 



BRG ERR 
3° 
3° 
3° 



Note from Table 2 that the reference point is at Station 1. 

Program outputs for Scenario 2 are given in List 2 on Page 
21. The two bearing estimates in the list are the reciprocal of 
the target's estimate of the bearing of Station 1. 
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Position Estimate 




Station 3 



Figure 1. A position estimate and a .33 confidence region. 

Table 1 gives the bearings of the indicated bearing 
lines and List 1 gives the bearing and the range of 
the position estimate from Station 3. 
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LIST 1 



OPTION? 

TPE=1 EST=2 $12=3 

PRB=4 CGH=5 

OPTION? 

1 

OBS BRS? 

TO 

$Tfi BR6? 

334 

STfl RN6? 

13598 
BP.6 ERR? 

4 

0B$ BRS? 

324 

ST ft BRS? 

58 

STfl RN6? 

11350 
BRS ERR? 

3 

OBS BR6? 

OPTION? 

2 

BR6= 359.51 
RN6= 19494.39 
END 

OPTION? 

3 

SIZE? 

PRB* 0.86 
SHJ* 1737.32 
DIR* 17.69 
SHI* 1232.96 
fl= 6T29444.91 
END 

OPTION? 

4 

PRB? 

.9 

SIZE* 2.15 
SHJ* 1864.12 
DIR* IT. 69 
SHI* 1322.94 
R= TT4T559.T7 
END 



OPTION? 

5 

OBS BR6? 

7 

ST ft BRS? 

0 

STR RH8? 

O 

BRS ERR? 

4 

OBS BRS? 

OPTION? 

2 

BRS* 9.15 
RNS= 19553. T8 
END 

OPTION? 

7 

SIZE? 

PRB* 0.86 
SHJ= 1712.95 
DIR* 12.48 
SHI* 1129.90 
fl= 6080418.49 
END 

OPTION? 

4 

PRB? 

,'5 

SIZE* 2.15 
$f!J= 1837.97 
DIR* 12.48 
SHI* 1212.36 
ft= 7900348. 38 
END 
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Figure 2. A position estimate and .71 confidence region. 

Table 2 gives the bearings of the indicated bearing 
lines and List 2 gives the reciprocal bearing and 
the range of the position estimate from Station 1. 
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List 2 



OPTION? 

Pi 


OPTION? 

c 


TPE=1 EST=2 SIZ=3 


.J 

OES ERG? 




326 


PRB=4 C0N=5 


STfi ERG? 




f>2 


OPTION? 

j 


STfi RHG? 

j TOCiu 


OBS ERG? 


1-jood 

ERG ERR? 


31 


T 


STfi ERG? 


OBS ERG? 


H 




STfi RHG? 


OPTION? 


8 


1 


ERG ERR? 

7 


BRS= 31.65 


OBS ERG? 


RNb= 13765.38 


352 




STfi ERG? 


END 


115 




STfi RNG? 


OPTION’ 


I110O 


3 


ERG ERR? 


SIZE? 


3 


2 


OBS ERG? 


PRB 5 9.86 


OPTION? 


Sfl.J- 2823.19 


2 

BR8= 31.99 


DIR= 172.86 


RHG- 14792.53 


SHI 5 967.90 


END 


fi 5 6151983.32 


OPTION? 


END 


T 




SIZE? 


OPTION? 


L 


4 


PRB= 8.36 


PRB? 




.9 


SHJ= 3696.56 


SIZE 5 2.15 


DIR- 15.39 


SHJ= 2178.35 


SMI = 1253.72 


DIR- 172.86 


fi= 14205070.9! 


SHI 5 1833.54 


END 


H 5 7882732.54 


OPTION? 


END 


4 




PRB? 

q 




SIZE 5 2.15 




SHJ= 3869.77 




DIR= 15.39 




SMI 5 1345.22 




fi 5 16354192.26 




END 
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IV. Program Listings 
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HP-4 1CV PROGRAM 



01*LBL "TPE" 

02 DEG 

03 CLRG 

04 .001 

05 STO 09 

0b 2b 

07 STO 10 
08+LBL -CON- 

09 FIX 0 

10 -OBS BRG?" 

11 PROMPT 

12 STO 11 

13 'STft BRG?‘ 

14 PROMPT 

15 STO 12 

16 ‘ST A RHG?" 

17 PROMPT 
13 STO 13 

19 -BRG ERR?" 

20 PROMPT 

21 STO 14 

22 1 

23 ST+ 03 

24 2 

25 RCL 03 

26 JOY? 

27 GTO 02 
23 1 

29 ST+ 10 

30 RCL 11 

31 STO IHD 10 

32 1 

33 ST+ 10 

34 RCL 12 

35 STO IND 10 

36 1 

37 ST+ 18 
33 RCL 13 

39 STO IHD 10 

40 1 

41 ST+ 10 

42 RCL 14 

43 STO IHD 10 

44 1 

45 RCL 03 

46 X=Y? 

47 GTO -COH- 

48 RCL 29 

49 RCL 28 

50 RCL 27 



101 


ST+ 


10 


102 


RCL 


IND I 


103 


STO 


12 


184 


1 




105 


ST+ 


19 


106 


RCL 


IND : 


107 


STO 


13 


108 


1 




109 


ST+ 


10 


110 


RCL 


IHD 


111 


STO 


14 


112 


XEQ 07 


113 


ISG 


@9 


114 


GTO 35 


115 


GTO 


■COH 


116»LBL 


02 


117 


XEQ 07 


113 GTO 


■COH 


119*LBL 


0? 


120 


RCL 01 


121 


RCL 


13 


122 


RCL 


12 


123 


SIN 




124 


* 




i *|C 

lij 


- 




126 


RCL 


02 


127 


RCL 


IT 

L J 


123 


RCL 


12 


129 


COS 

* 




i -*FJ 

131 


•> 




132 


R-P 




133 


STO 


16 


134 


v/ y 


i 


135 


STO 


15 


136 


RCL 


11 


137 


RCL 


4 C 

1 J 


133 


- 




139 


180 




140 


X<>' 


f 


141 


V 


Y? 


142 


GTO 


11 


143 


360 




144 


- 




145»LBL 


11 


146 


RCL 


14 



147 / 

10 148 STO 17 

149 RCL 16 

150 RCL 14 



51 - 

52 SIH 

53 * 

54 STO 15 

55 RCL 33 

56 RCL 32 

57 RCL 31 

58 - 

59 SIH 

68 * 

61 STO 16 

62 RCL 31 

63 RCL 27 

64 - 

65 SIH 

66 STO 17 

67 X=0? 

63 GTO 08 

69 RCL 15 

78 RCL 31 

71 SIH 

72 * 

73 RCL 16 

74 RCL 27 

75 SIH 

76 * 

77 - 

78 RCL 17 

79 / 

30 STO 01 

81 RCL 15 

32 RCL 31 

33 COS 

OA ± 

35 RCL 16 

36 RCL 27 

37 COS 

83 * 

89 - 

90 RCL 17 

91 / 

92 STO 82 

93 26 

94 STO 10 

95*LBL 05 

96 1 

97 ST+ 18 

93 RCL IHD 

99 STO 11 

100 1 



23 



151 * 

152 PI 

153 * 

154 130 

155 / 

156 STO 18 

157 RCL 15 
153 COS 
159 RCL 13 
168 / 

161 STO 19 

162 RCL 15 

163 SIN 

164 RCL 18 

165 / 

166 STO 20 

167 RCL 19 
163 Xt2 

169 ST+ 03 

170 RCL 19 

171 RCL 20 

172 * 

173 ST+ 04 

174 RCL 20 

175 Xt2 

176 ST+ 05 

177 RCL 17 

178 RCL 19 

179 * 

180 ST+ 06 

131 RCL 17 

132 RCL 20 

133 * 

184 ST+ 07 
135 RTH 
136»LBL -EST- 
137 FIX 2 
183 RCL 04 

189 Xt2 

190 RCL 83 

191 RCL 85 

192 * 

193 - 

194 STO 20 



195 


X=0? 


196 


GTO 


08 


197 


RCL 


01 


198 


RCL 


04 


199 


RCL 


07 


200 


* 





251 


COS 




252 


Xt2 




253 


RCL 


05 


254 


* 




255 


RCL 


23 


256 


1 




257 


P-R 




253 


* 




259 


2 




268 


* 




261 


RCL 


84 


262 


* 




263 


STO 


26 


264 


- 




265 


RCL 


23 


266 


SIN 




267 


Xt2 




268 


RCL 


03 


269 


* 




278 


+ 




271 


RCL 


20 


272 


CHS 




273 


/ 




274 


STO 


24 


275 


RCL 


23 


276 


SIN 




277 


Xt2 




273 


RCL 


05 


279 


* 




280 


RCL 


26 


231 


+ 




282 


RCL 


i-%1 


283 


COS 




284 


Xt2 




285 


RCL 


03 


236 


* 




237 


+ 




233 


RCL 


20 


289 


CHS 




298 


/ 




291 


STO 


25 


292 


RCL 


24 


293 


XOY? 


294 


GTO 


01 


295 


STO 


25 


296 


XOY 


297 


STO 


24 


293 


90 




299 


ST+ 


23 


300»LBL 


01 



261 RCL 85 

202 RCL 06 

293 * 

294 - 

205 RCL 20 

206 / 

207 + 

203 STO 26 

209 RCL 02 

210 RCL 03 

211 RCL 07 

212 * 

213 RCL 04 

214 RCL 06 

215 * 

216 - 

217 RCL 20 

218 / 

219 + 

220 RCL 26 

221 XOY 

222 R-P 

223 STO 21 

224 XOY 

225 X>0? 

226 GTO 06 

227 360 

223 + 

229*LBL 06 

230 STO 22 

231 RCL 04 

232 RCL 03 

233 RCL 85 

234 X*Y? 

235 GTO 03 

236 RCL 04 

237 SIGN 

233 45 

239 * 

240 GTO 84 

24HLBL 03 

242 - 

243 / 

244 2 

245 * 

246 AT AH 

247 2 

248 / 

249*LBL 04 

250 STO 23 



24 



381 -BRG=' 

392 flRCL 22 
363 ft V I EM 

304 STOP 

305 -RNG=* 

30b flRCL 21 

307 flVIEH 

308 STOP 
389 -EHD- 
310 GTO 80 

31 1 ♦LBL -SI2- 

312 -SIZE- 

313 PROMPT 

314 STO 11 

315 XT2 
31b 2 

317 / 

318 CHS 

319 EtX 

320 1 

321 XOY 

322 - 

323 STO 80 

324 *PRB=- 

325 flRCL 08 

326 flVIEH 

327 STOP 

328 GTO 09 
329»LBL -PRB- 

330 - PRB? • 

331 PROMPT 

332 CHS 

333 1 

334 + 

335 LH 

336 2 

337 * 

333 CHS 
339 SORT 
348 STO 11 

341 -SIZE=- 

342 flRCL 11 

343 flVIEH 

344 STOP 
3454LBL 09 

346 RCL 23 

347 STO 12 

348 RCL 11 

349 RCL 25 
358 SORT 



351 * 

352 RCL 11 

353 RCL 24 

354 SORT 

355 * 

356 XOY 

357 -SMJ=- 

358 flRCL X 

359 flVIEH 

360 STOP 

361 RCL 23 

362 X>8? 

363 GTO 18 

364 130 

365 + 

366*LBL 18 
367 " DRC= ' 
363 flRCL X 
369 flVIEH 
378 STOP 

371 -SMI=- 

372 RBH 

373 flRCL Y 

374 flVIEH 

375 STOP 

376 PI 



379 -fl=- 
389 flRCL X 
381*LBL 80 
382 flVIEH 
333 STOP 
384 -END" 

335 GTO 08 
386*LBL 83 

387 flVIEH 

388 STOP 

389 -NO SOL- 
398 GTO 08 
391 EHB 
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SHARP PC-1500 PROGRAM 



10: "A": CLEAR :Din 
A (7) : DEGREE 
15: "C": INPUT "OBS 
BRG? " ; P : 
INPUT "STA BRG 
? " J Q: INPUT "S 
TA RNG? "JR: 
INPUT "BRG ERR 
? "JO: PAUSE " 

II 

20: IF I =2G0T0 50 
25: 1=1+1 : A( I- I )=P 
: A< I +1 )=Q: A( I + 
3>=R: A( l+5)=0: 
IF I = 1 GOTO 15 
30: X=A(4)YSIN (A( 
2)-A(0) ) : Y=A(5 
)*SIN (A(3)-A( 
1 ) ) : Z=S I N CA ( 1 
)-A<0) ) : IF Z=0 
GOTO 150 

35: U=(XYS1N A(1 )- 
YYSIN A(0))/Z: 
U=(XYCOS A ( 1 ) - 
YYCOS A(0))/Z 
40: FOR n=0TO 1 
45: P=A(f 1 ) : Q=A<n+2 
) : R=A <(1+4 ) : 0=A 
(H+6) : GOSUB 17 
0 : NEXT M: GOTO 
15 

50: GOSUB 170: GOTO 
15 

60: "Z": PAUSE " ": 
F=(BYB-AYC) : IF 
F=0GOTO 150 
65: X=U+ (BYE-CYD )/ 
F: Y=U+ ( AYE-BYD 
)/F: GOSUB 200 
70: T =SGN BY45: IF 
A=CGOTO 80 
75: T=. 5YATN (2*B/ 

< A-C ) ) 

80:G=(C*COS TYCOS 
T-2YBYC0S TY 
SIN T+A*SIN T Y 
SIN T )/-F : G=TG 
85: H=(CYSIN TYS1N 
T+2YBYC0S TY 
SIN T+AYCOS Tt 
COS T)/-F:H=TH 
: IF H>=GGOTO 3 
5 

30: Z=H: H=G: G=Z: T= 
T+90 



95: PRINT "BRG="JJ 
.•PRINT "RNG=" J 
K : GOTO 145 

100: "H": PRINT " TP 
E=A EST=Z SIZ= 
S PRB=X" 

105: PRINT "CON=C" : 
GOTO 100 

120: "S": INPUT "SJZ 
E? " J S: 0=1-EXP 
(-SYS/2) 

1 25 : PR I NT USING "tt 
tt. tttt"; "PRB="J0 
: US I NG : GOTO 1 
35 

130: "X" : INPUT "PRB 
? " J 0: S=T (-2Y 
LN (1-0)): 
PRINT "S1ZE="J 
S 

135: X=SYG: Y=SYH: N= 
T.'PRINT "Snj = " 
J Y: IF N<0LET N 
=N+ 1 80 

1 40 : PR I N T "DJR="JN 
: PR I NT "Sni = "J 
X : PR I NT "A="jn 
YXYY 

1 45 : PR I NT "END": 
GOTO 145 

150: PRINT "NO SOL" 

: GOTO 150 

170: X=U-RYSIN Q:Y= 
U-RYCOS Q: 
GOSUB 200 

175: U=(P-J) : L=KYOY 
JI/180: G=COS J/ 
L: H=S1N J/L: IF 
Ui >= 1 80LE T W=W- 
360: GOTO 185 

180: IF U<=- 1 80LET 
W=W+360 

185: UJ=U/0: A=GYG+A: 
B=GYH+B: C=HYH+ 
C : D=WYG+D : E=WY 
H+E: RETURN 

200: K=T (XYX+YYY ) : 

IF K=0LET J=0: 
RETURN 

205: J=ACS (Y/K) : IF 
ASN (X/KX0LET 
J=360- J 

210: RETURN 
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SHARP PC-1211 PROGRAM 



10: "A" : CLEAR : 
DEGREE 

15: "C": INPUT "0 
BS BRG? "iP: 
INPUT "STA B 
Rij? ">Q: 
INPUT “STA R 
HG? "iR: 
INPUT "BRG E 
FR? "50: IF I 
=260 TO 130 

75: 1=1+ 1!A< 1+26 
>=P: A< I+28)= 
Q: A< 1+30 >=R: 
ACI+32>=0: IF 
1=1 GOTO 15 

35:X=AC31>*SIN 
CAC280-A C27) 
>: Y=AC32 >* 
SIN CAC301-A 
C28)>:Z=S1N 
CAC28)-AC27> 

90: IF Z =060 TO 3 
25 

95: U=<X*SIN A <2 
8>-Y*<3IM AC 2 
V=CX* 
COS AC 23 >-Y+ 
COS AC 27 IO/2 
105: FOR 1=1 TO 2: 
F-ACI+26>:Q= 
AC 1+28) : R=A< 
1+30 > : 0=AC 1+ 
32 ) : GOSU B 40 
0 

110: NEXT I: GOTO 
15 

130: GOSU B 400 
135: GOTO 15 
140: "Z" : F=CB*B-A 
*€> : IF F=0 
GOTO 325 

150: X=U+ CB+E--C+B 
>/F : V-V+ CA+E 
-B+DX-'F: 
GOSUB 500 
160: T=SGH B*45: 

IF ACTHEN 1 
70 

165: T=. 5+ ; ATN <2* 
£■•••'( A -C> > 
170:G=<C*C0S T+ 
COS T-2+-E+ 
COS T*SIN T+ 
A+SIN T+SIN 
T)/-F:6=-TG 
175: H=0>SIN T* 
SIN T+2*B* 
COS T+SIN T+ 
A+COS T+COS 
T>/-F:H=TH: 

IF H>=GGOTO 
185 



130: Z=H: H=G = 6=Z: 
T=T+90 

135: PRINT "BRO= " 
5 J: PRINT "PH 
G=" j K: GOTO 3 
20 

200: "H": PRINT 11 T 
PE=A EST=Z S 
IZ=S PRB=X" 

205: PRINT "CON=C 
"•.GOTO 200 

225: "S": INPUT 
IZE? " i S: 0=1 
-EXP C-S+S/2 

230: PRINT USING 
tttt" i "F'RB 
="50: USING : 
GOTO 300 

235: "X": INPUT "P 
RB? " 5 0: S=T C 
-2+LH C 1 — 0 ) ) 

: PR I NT "SIZE 
= ";s 

300: x=S*G: Y-S+H: 
N=T: PRINT "S 
PI J=" ; Y: IF N 
<=OLET N=N+1 
80 

315: PRINT "DRC=" 
SN: PRINT "SM 
I = " i Z=Tf+X+ 

V: PRINT "A= H 

• 

? < 2 . 

320: PRINT "END": 
GOTO 320 

325: PRINT "NO SO 
L" •' GOTO 325 

400: X=IJ-R+SIN Q: 
V'=y-R'+ : CuS Q: 
GOSUB 500:W= 
C P-..T > : L= K+0+ 
jr/130: G-COS 
J.-'L : H=S I N Jy 
L 

415: IF U >= 1 8 OL.ET 
W=M~36ij: GOTO 
440 

425: IF NO- 180 
LET W=W+360 

440: H=N 0: A=G+G+ 
A: B=G+H+B: C= 
H+H+C: H=H+G+ 
D:E=H*H+E: 
RETURN 

500: k=T CX+X+Y+Y) 

: IF K=OL ET J 
=0: RETURN 

510: M= AS N C X .-'K > : 
J=ACS CY -'I-O: 
IF M COLEY J= 
360- J 

515: RETURN 
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CASIO FX-702P PROGRAM 



5 INF "OPTION"! J: 
SET F2i MODE 4 
18 IF M THEN 175 

15 If J=2 THEN 180 

28 IF J=3 THEN 145 

25 IF J=4 THEN 150 

30 IF J=5 THEN 45 

48 VflC ’.SET F2 
45 INF ”068 BR6',P 
i INF "STfi BR6”, 
Q: INF ”$Tfi RN6” 
,R 

58 INF ”BR6 ERR". 0 
:IF 1=2 THEN 95 
55 I=I+i:fl(I-l)=P: 
fi(I+l)=Q:fi(I+3) 
=R:fi(I+5)=0 
68 IF 1=1 THEN 45 
65 X=fi(4)*SIN (fi(2 
H(0)):y»fl(5)* 
SIN <fl<3H(l» 
79 Z=SIH (fi( 1 )-R(9 
)): IF 2=6 THEN 
170 

75 IJ=(X*$IN HUM 
♦SIN fl(0))/Z:V= 
<X*CGS R(1H*C 
OS fi(8))/Z 
88 FOR H=8 TO l:P= 
R(H):6=fl(N+2) 

35 R=fi(N+4):ij=fl(N+ 
6):SSB 288: NEXT 
M: SOTO 45 

99 R=fi(M+4):0=fi(M+ 
6)sfiSB 290! NEXT 
H! SOTO 45 

95 SSB 208: SOTO 45 
188 F=(B*6-fi*C):lF 
F=6 THEN 170 
185 X=U+(B*E-C*D)/F 
:Y=V+(fi*E-6*D>/ 
F: RPC Y,X 

118 T=S6N 6*45: IF H 
=C THEN 120 
115 T=.5*flTH (2*8/( 
fi-O) 

120 S=(C*COS T*COS 
T-2*8*C0S T*SIN 
T+fl*SIN T*SIH 
T)/-F 



125 H=(C*SIN T*SIH 
T+2*0*COS T*SIN 
T+fl+C-OS T*COS 
T)/-F 

138 6=S8R S;H=8SR H 
••IF i3>H;Z=H:H=S 
:8=Z:T=T+90 

135 IF y<0;Y=Y+368 

149 FRT ”8RS=”;Y:PR 
T "RH6=*;X:G0T0 

165 

145 INF ’’SIZE” .8:0= 
1-EXP (-8*8/2): 
FRT •’•PRB=”;0:GQ 
TO 155 

159 INP ”PRBMj:S=S 
OR (-2+LH (1-0) 
):RRT ”SIZE=-”;S 

155 X-S*6:Y-S*H:PRT 
”SHJ=”iY! N=T: I 
F N(8;N=f!tl80 

166 PRT *DIR="?N:PR 
T "SM=";X:PRT 
”fi-”‘!I*X*T 

165 FRT "END”: SOTO 
165 

179 PRT "HO SOL”: 80 
TO 170 

175 PRT ”TFE=1 EST= 
2 817 = 5 ” 

150 PRT *PRB=4 CON= 
5”: SOTO 175 

290 X=U-R*SIH Q:\-V 
-R*COS Q : RPC Y> 
X 

295 N=<F-Y):L=X*0*i 
/ i 80 : S=COS Y/l: 
H=SIN Y/L 

218 IF H±1885H=H-36 
8 : GOTO 226 

215 IF Hi-189; H=H+3 
68 

228 H=W/0:fi=e*6+R:B 
=S*H+B:C=H*H+C: 
D=M*6+D:E=H*H+E 
iRET 
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TI-59 PROGRAM 



000 


47 


CMS 


050 


10 


E 3 


1 00 


42 


STO 


001 


6 0 


DEG 


051 


18 


C 3 


101 






002 


02 


il 


052 


p / 


IFF 


102 


43 


RCL 


003 


42 


STO 


053 


03 


0 3 


103 


14 


14 


004 


09 


03 


054 


01 


01 


104 


39 


CDS 


005 


69 


OP 


055 


28 


2 o 


105 


65 


x 


006 


«- 

iL O 


2 y 


0 5 6 


13 


D 3 


106 


43 


RCL 


007 


76 


LBL 


057 


75 


- 


107 


18 


18 


008 


16 


R 5 


058 


13 


D 3 


108 


75 


- 


003 


03 


Q 


053 


35 


= 


103 


43 


RCL 


0 1 0 


42 


STO 


060 


34 


+ /- 


1 1 0 


10 


10 


Oil 


0 0 


00 


061 


•j o 


S I H 


1 1 1 


39 


COS 


012 


42 


STD 


062 


65 




1 12 


65 




0 1 3 


01 


01 


063 


13 


D 3 


1 1 3 


43 


RCL 


0 1 4 


9 2 


RTH 


064 


35 


= 


114 


13 


13 


015 


•• 

( h: 


LBL 


065 


48 


EXC 


115 


35 


= 


0 1 6 


1 8 


C 1 


0 b 6 


13 


13 


116 


55 


-i- 


0 1 7 


69 


OP 


Ob? 




I HV 


117 


43 


RCL 


0 1 8 


20 


20 


068 


97 


DSZ 


118 


20 


20 


0 ! 3 




ST* 


069 


09 


03 


1 13 


“j Q 


S I H 


020 


0 0 


00 


070 


0 0 


00 


1 20 


35 


= 


021 


Q *? 

•• L- 


RTH 


07 i 


f' b 


7 6 


121 


42 


STD 


022 


7 6 


LBL 


U i LI 


01 


1 


1 22 


39 


33 


U ti o 


13 


D 3 


073 


63 


OP 


123 


86 


STF 


024 


63 


OP 


0/4 


21 


21 


124 


0 3 


03 


025 


21 


21 


075 


32 


RTH 


125 


02 




026 


73 


PC* 


076 


42 


STO 


126 


42 


STD 


ij v 7 


01 


01 


077 


18 


18 


127 


03 


03 


028 


3 2 


RTH 


078 


65 


X 


128 


16 


fi 3 


023 


7 6 


LBL 


0/3 


43 


RCL 


123 


13 


D 3 


030 


10 


E 3 


080 


14 


14 


130 


10 


E 3 


031 


65 


V* 


081 


42 


STO 


131 


42 


STD 


032 


3 3 


Vi 


082 


20 


20 


132 


18 


18 


033 


55 


-r 


083 


*“ s I~2 


S I H 


133 


13 


D 3 


034 


0 1 


i 


084 


75 


- 


134 


.3 ^ 


X i T 


035 


08 




085 


43 


RCL 


135 


13 


D 3 


036 


00 


0 


086 


10 


10 


136 




X i T 


n 3 7 


35 


= 


n i-i 7 


.•* .«* 


I HV 


137 


: 


P/R 


038 


32 


RTH 


088 


44 


SUM 


138 


75 


- 


033 


7 6 


LBL 


ij y 9 


20 


20 


133 


43 


RCL 


040 


l i 


H 


09 0 


. -| i"« 


S I H 


140 


•J i“: 




041 


81 


RST 


09 i 


65 




141 


35 


— 


042 


7 6 


LBL 


092 


43 


RCL 


142 


34 


4- / - 


043 


12 


B 


0 93 


13 


13 


143 




X I 


044 


i y 


C 3 


094. 


35 


— 


144 


75 


- 


045 


31 


R/S 


095 


55 


-r 


145 


43 


RCL 


046 


i y 


c 3 


096 


43 


RCL 


146 


3 9 


39 


047 


31 


R/S 


097 


20 


20 


147 


35 


= 


048 


18 


C 5 


098 


i": 

•Z l *- J 


S I H 


148 


34 


+ / - 


043 


31 


R/S 


099 


35 


= 


143 




K 1 1 



29 



150 


# ~t 


IHV 


151 


I s f' 


P/R 


152 


42 


S TO 


153 


19 


19 


154 


10 


E * 


155 


75 


- 


15b 


43 


RCL 


157 


18 


18 


153 


95 


— 


159 


94 


+ / - 


160 


3 


X i T 


161 


69 


OP 


1 62 


21 


21 


i 6 3 


64 


PD* 


164 


01 


01 


165 


42 


STO 


166 


18 


18 


167 


89 


tf 


1 6 y 


O ul 


X 1 T 


169 


"i "i 


INV 


170 


77 


GE 


171 


01 


01 


172 


~y i” i 
j* «“i 


78 


173 


-? nr 
i J 




174 


•"i •”« 


x i i 


175 


65 


x 


176 


02 


V 


177 


95 


= 


178 


49 


PRD 


179 


18 


18 


180 


73 


EC* 


1 8 1 


01 


01 


182 


35 


1 /X 


183 


42 


STO 


184 


•j 7 


27 


185 


42 


STO 


186 


26 


26 


187 




XI T 


188 


43 


RCL 


189 


19 


19 


190 


37 


P/R 


191 


42 


STO 


192 


' J »_ s 
«■ i_i 


28 


193 


4y 


PRD 


194 


, —j 


• J ~7 


195 


j 


E 


196 


44 


SUM 


1 9 !* 


ill • J 


23 


198 


43 


RCL 


199 


18 


18 



200 
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Appendix 1. The Estimation Procedure 

In this development the assumptions stated in Section I are 
required conditions. A rectangular coordinate system is used 
with the positive y-axis directed north, the positive x-axis 
directed east and the origin at the reference point and all angles 
are in radians. Figure 3 shows a station located with respect to 
the coordinate system. The bearing line of length r goes to the 
object's unknown position, the bearing line of length r goes to 
an initial estimate of the object's position and the third bear- 
ing line corresponds to an observed bearing. 

To determine the coordinates for a final estimate, consider 
the arc segments u = r ( 0 — cj) ) where 0-<J> is the bearing error and 
v = r(cj)-£) and w = r(0-8) that are defined by the three bearing 
lines and the circle of radius r that is centered on the station 
and goes through the initial estimate. The geometry is shown in 
Figure 3. Note that u can be defined by u = w - v. In this 
expression, w = r(0-B) is known, and v can be determined in terms 
of x and y the unknown coordinates of the target. With the refer- 
ence point at the initial estimte, to first order, 
v = x cos 8 - y sin 8 and u = w - x cos 8 + y sin 8. Since this 
approximation applies to all stations, its use suggests that, 
for each station, r z r which is equivalent to having the initial 
estimate relatively close to the target's position. 

Since, for each station i, an observed bearing 0^ is the 

value of a normal random variable 0^ with mean <t>^, the coordinate 

u. = r.(0.-<J>.) is the value of a normal random variable U. with 

1111 1 

mean zero. In addition, since the 0 i are independent, the tb are 
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Figure 




. The coordinate geometry. The coordinates of the initial 
estimate are (x*,y*). In the development, the reference 
point is at the initial estimate. 
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also independent. And, consequently, the joint distribution of 
the It is determined. 

To estimate x and y, maximum likelihood estimates are used 
here. The likelihood for a sample 0^, 0 2 ,..., 0^ from n stations 
is 

^1 I n o 9 

L(0, ,0.,. ../0 ) = II exp - Z (0.-cf>.) /e . 

12 n 1 /SW e. 2 1 i 1 1 



and the likelihood for a cooresponding sample u^, u 2 ,..., u^ is 



L(u. ,u 0 , u ) = 

1 z n 



n i 

n — - — 

1 /2tt a . 

l 



1 n 2 2 
exp - j Z u I/O' 



with a. = r.e. where e. is the standard deviation of 0.. 
1111 i 

By definition, the maximum likelihood estimates for x and y 
are the estimates which make L(u^, u 2 ,...,u n ) a maximum. In 
this case, making L (u^ , u 2 , . . . , u n ) a maximum is equivalent to 



n 



2 , 2 



making Z(u. /a. ) a minimum. So, to find the maximum likelihood 



l ' l 



estimates x and y, solve the following two equations for x and y: 



3( * nL) .| =0 and 3( ! n L) | =0 



3x 



3y 



x = x 



= y 



With u^ = w^ - x cos 3^ + y sin 3^, the equations can be written 
as follows: 



n 



E (w^ - x cos 3 i + y sin 3^) (cos 3 ^)/ck = 0 



and 



^ /S /N , 2 « 

E (w^ - x cos 3^ + y sin 3^) (sin 3j_)/cK = 0 
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In terms of the following quantities: 



A = E (cos^ , B = £(sin 3^ cos B^)/o? , 

C = £(sin 2 3 i )/a 2 , D = £ (w^ cos 3 i )/a? , 

E = £ (w i sin 3 i ) /a? , 

the equations become: 

Ax - By = D 

Bx - Cy = E 



The solutions are: 

(1) x = (BE - CD)/(B 2 - AC) 
and 

(2) y = (AE - BD)/(B 2 - AC) 

A confidence region can be constructed about an estimated 
position. In order to indicate how this can be done, a proba- 
bility region about the true position will be considered first. 

Both x and y are values of random variables. If a new set 
of bearings 0^, & 2 ' • • • / 0 n is observed (for the same initial 
estimate and a fixed target) , in general, a new pair of values 
x and y will be obtained. 
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If 



X and Y represent these random variables, 



/\ 

X 



1 

(b 2 -ac) 



£ (W- /°f) (B sin 3; 
1 1 1 1 



C cos 3^) 



Y 



1 

(b 2 -ac) 



2 (W./O 2 ) (A sin 3, 
1 1 1 1 



B cos 3^) 



where W. = r.(0. - 3-). 

l i i i' 

A /N 

Since X and Y are a linear combination of the n normal random 
variables W-^, W^, or equivalently of the n normal random 

variables 0^, © 2 , © n , they have a joint normal distribution. 

Since E(W i ) = r ± ( 4> ± - 3 i ) , if 3 i = 4> i for i = 1, 2, n, that 

is, if the initial estimate of the target's position is at the 
target's position, E(W 1 ) = 0 for i = 1, 2, ..., n. In this case 
E(X) = 0 and E(Y) = 0 and the joint normal distribution is centered 
on the object's position. To the degree of the approximations 
that have been made, this is also true if the initial estimate is 
not at the target's position. 

A region of minimum area for a given probability of contain- 
ment of an estimated position can be determined. The region is 
bounded by an ellipse which is centered on the object's position 
and whose axes lie along the axes of an x ' y ' -coordinate system 
that is obtained by rotating the xy-coordinate system that is 
centered on the object's position through an angle y. In this 

system, a-', ~ , is 0, that is x' and y' are independent normal ran- 
x y 

dom variables. The two coordinate systems are illustrated in 
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Figure 4. The coordinates of a point in the two systems are 
related by 

x' = x cos Y ~ y sin Y 

y' = x sin Y + y cos Y 

These relations imply: 

2 2 2 2 2 

(3) cu = cos Y ~ 2a cos y sin y + sin y , 

x ' x xy y 

2 2 2 2 2 

(4) = a^ sin y + 2a AA cos y sin y + cos y 

y ' x xy y 

and 

2 2 2 2 2 

(5) ^ = ( a ^ - a^) sin y cos y + a^ (cos y - sin y) 

x ' y ' x y xy 



where y, the angle of rotation of the coordinate axes, is posi- 
tive in the clockwise direction. And = 0 implies 

x ' y ' 



tan 2 y 



2a~^ 

xy 




y x 

With the initial estimate of the target's position at the 

2 

target's position E(W^) = 0 and therefore, Var(W^) = a^ for 
i = 1, 2, ..., n. In this case 



o ^ E (1/a?) (B sin (3. - C cos 6-) 2 , 

x (B -AC)^ 1 1 1 1 
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North 




East 

Figure 4. An elliptical confidence region and the primed 
coordinate system in which the covariance 
is zero. The center of the ellipse and the 
origin of the coordinate systems are at the 
target's estimated position. The estimated bear- 
ing a and estimated range p are indicated for a 

reference position. 
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1 



2 

a = 



2 2 
( B -AC) 



n 

Z 

1 



(1/a?) (A sin g^ - B cos g^) 2 



and 



2 

= — 5 2 (l/a.)(B sin g. - C cos g.)(A sin g. - B cos g . ) . 

xy (B -AC) 1 1 1 11 1 



Using the definition for A, B and C, the above become 



( 6 ) 

(7) 



and 



„2 _ C 

0 ~ 2 ’ 

x (AC-B ) 



/N 2 ' 

y (AC-B ) 



( 8 ) 




(AC-B 2 ) 



So tan 2 = 2B/(A-C) for g^ = <f>^, i = 1, 2, . n. 

With the target's position known and, consequently, <|k known 

2 2 

for i = 1, 2, ..., n, the above expressions for Og , , a ~~ and y 

can be used, since the initial estimate of the target's position 

can be taken as the target's position. 

With values for , o ~ and y, values for a~ and o* 

x y xy x y ' 

can be found by using equations (3) and (4). The probability 
that an estimated position will be within an ellipse of semiaxes 
ko£, and kc~, which is centered on the target's position is 
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1 -exp ( -k /2) . This result can be found by integrating the bi- 
variate normal density over the ellipse. And the area of the 
2 

ellipse is k' a^,a~. 

x' y ’ 

Given estimates x and y found by using Equations (1) and (2) , 
the ellipse with semi-axes ko~, and ko~, in a coordinate system 



A A 

that is centered on the point (x,y) and has been rotated through 

2 

an angle y is a 1 - exp(-k /2) confidence region. This follows, 
since, to the degree of the approximations involved, the bivar- 
iate normal distribution of X and Y is centered on the target's 

2 2 

position. The confidence ellipse is defined if o~ , and 

c x y xy 

can be found, that is if the elements of the covariance matrix 
can be found. To the degree of the approximations involved, 
this can be done as follows: First, assume the initial estimate 

of the target's position is at the target's position. Then, 

2 2 . . 
values for and y can be determined by using Equations 

x > y i xy ^ 

(6), (7) and (8). These values can then be used to determine 

2 2 

a~, , a-' , and a~,~, by using Equations (3), (4) and (5). Now, 

x y x y 

with a value for k, a confidence region can be constructed. To 
the degree of the approximations involved, the shape of the con- 
fidence region is independent of both the target's position and 
of the initial estimate of the target's position. 

For the case where bearings are taken from the target on two 
or more stations, 0 i is the reciprocal of the bearing taken from 
the target . 

A discussion of the theory of bearings only position 
estimation procedures for situations similar to the one considered 
here is given in Reference 1. Reference 2 gives an essentially 
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equivalent bearings only procedure. It also gives a range only 
procedure, a range and bearing procedure and HP-9830A programs 
with which to 'implement the procedures. Using the fix determined 
by two lines of bearing as the initial estimate was suggested by 
this reference. 

The following equations are evaluated in the program to deter- 
mine the coordinates of the initial estimate: 

x* sin (82 - 9^) = [p-^ sin (a-^-0-^) ] sin 02 

- [p 2 sin (a 2 - 02 )] sin 0-j_ 

and 



y* sin ( 02 ~ 0 -^) = [p^ sin (a^-0^) ] cos ©2 

- [p 2 sin (a 2 ~ 02 ) 1 cos 6^ • 

Reference 3 describes a TI-59 program that is based on an equi- 
valent procedure. The program allows a user to either input 
the coordinates of the initial estimate or determine them in 
the manner described here. 
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Appendix 2. HP-41CV Program Labels 

The global labels in the HP-41CV program that are assigned 
to the keys Z+, 1/x, /x, LOG and LN give these keys a mnemonic 
character. For example, with the calculator in USER mode, if 
LOG is pressed and held, SIZ will be displayed and then after 
a delay, NULL. 

If a global label is replaced by a local label, this mnemonic 
character will be lost. However, automatic key assignment will 
be gained if the key corresponding to the local label has not 
been previously assigned. If this is the case and if there is 
a second program in program memory that uses the same local label, 
then it will be automatically assigned when the calculator is 
positioned in program memory at that program. 

Global labels in the HP-41CV program can be replaced with 
local labels. If this is done as described below, the user in- 
structions will still be applicable. First, make the following 
replacements : 



Line 


329: 


LBL 


" PRB " 


-> 


LBL 


E. 


Line 


311: LBL "SIZ" 


-> 


LBL D. 


Line 


186: 


LBL 


"EST" 


-> 


LBL 


B . 


Line 


008: LBL "CON" 


-> 


LBL C. 


Line 


047: 


GTO 


"CON" 


-> 


GTO 


C. 










Next , 


after 


line 


001 : 


LBL 


"TPE" 


, insert LBL A 


SO 


that it 



becomes line 002. Finally, if there is a key assignment for a 
key in the following list: £ + , 1/x, /x, LOG and LN, remove it. 

Now, with the calculator in USER mode, if LOG is pressed and held, 
XEQ D will be displayed and then, after a delay, NULL. The key's 
mnemonic character has been lost, but if there is a second pro- 
gram in program memory with the local label D, LOG will be auto- 
matically assigned to D when the calculator is positioned in pro- 
gram memory at that program. 
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